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The implementation of the linear covariant gauge on the lattice faces a conceptual problem: using 
the standard compact discretization, the gluon field is bounded, while the four-divergence of the 
gluon field satisfies a Gaussian distribution, i.e. it is unbounded. This can give rise to convergence 
problems when a numerical implementation is attempted. In order to overcome this problem, one 
can use different discretizations for the gluon field or consider an SU(iV c ) group with sufficiently 
large N c - One can also consider small values of the gauge parameter £, and study numerically the 
limiting case of B, — > 0, i.e. the Landau gauge. These different approaches will be discussed here. 
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1. Why study the linear covariant gauge? 

Extensive numerical simulations on very large lattices ||, ||, ||, \j§ have shown that, in 
three and in four space-time dimensions, the Landau gluon propagator displays a massive solution 
at small momenta and the Landau ghost propagator is essentially free in the same limit. (For a re- 
cent review and various comments, see respectively ^ and [§].) These results are not in agreement 
with the original Gribov-Zwanziger confinement scenario [[(], pit] but they can be explained in the 



so-called refined Gribov-Zwanziger framework [12, 13]. 



A natural extension of these works in Landau gauge would be to consider the linear covariant 
gauge, which has the Landau gauge as a limiting case. However, until recently, the numerical 
gauge fixing for the linear covariant gauge [ p], p3| , 16] was not satisfactory. In Reference [17] 



we introduced a new implementation of the linear covariant gauge on the lattice that solves most 
problems encountered in earlier implementations. The final goal is to evaluate Green's functions 
numerically in the linear covariant gauge and, in particular, in the Feynman gauge. This could 
allow a nonperturbative evaluation of the gauge-invariant off-shell Green's functions of the pinch 



technique [18, 1 



2. The Linear Covariant Gauge on the Lattice 

In the continuum, the linear covariant gauge is obtained by imposing the gauge condition 

d^A h l {x)=A b (x), (2.1) 
where the real- valued functions A b (x) are generated using a Gaussian distribution 



exp 



1 

'2f 



d d xY, W(x) 



(2.2) 



with width y^. The limiting case t, — > 0, which implies A b {x) = 0, yields the Landau gauge. 
On the lattice, one can fix the Landau gauge by minimizing the functional 

<?lg[U s ] = -Tr£s(x)£/ M (x)s t (x + ^) (2.3) 

with respect to the gauge transformations {g(x)}. Here, U^(x) are (fixed) link variables and g{x) 
are site variables, both belonging to the S\J(N C ) group. The sum is taken over all lattice sites x and 
directions }X, and Tr indicates trace in color space. By considering a one-parameter subgroup 1 



exp 



izy b (x)X t 



(2.4) 



of the gauge transformation {g(x)}, one can verify that the stationarity condition for the functional 
<%lg [U 8 ] implies the (lattice) gauge condition 



£4(x)-Aj(x-^)=0, 
I 1 



(2.5) 



'Here we indicate with X b a basis for the SU(Af c ) Lie algebra and with '/'(x) any real-valued function. 
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where 



[li)- 1 [U^x)-Ul(x)] i 



(2.6) 

I traceless 

is the usual lattice discretization for the gluon field. Also, from the second variation of S'lgW 8 ] 
one can obtain a discretized version of the Faddeev-Popov operator 



ab 



-T) ab rl 



(2.7) 



where D ab is the covariant derivative. Clearly, for the gauge-fixed configurations, i.e. for all local 
minima of the functional <§lg [U g ] , this operator is positive-definite. This set of local minima defines 
the first Gribov region £2 [JToJ, [TTJ] . 

In Reference [17] we have introduced the minimizing functional 2 



(2.8) 



which is a natural extension of the Landau functional (2.3). Here 9^ indicates real part. One 
should stress that, in the numerical minimization, the link variables (x) are gauge-transformed 
to g(x)Ufi(x)g^ (x + e^, while the A b (x) functions do not get modified. It is easy to verify, using 
again a one-parameter subgroup g(x, t), that this functional leads to the lattice linear covariant 
gauge condition 



£A*(jc)-A*(x-« m )=A*(x). 



Note that the above relation and periodic boundary conditions yield 

2>*(*) = o. 



(2.9) 



(2.10) 



This equality must be enforced numerically, within machine precision, when the functions A b (x) 
are generated using the Gaussian distribution ( J2.2[ ). Also note that the second variation (with respect 
to the parameter t) of the term ig(x)A(x) is purely imaginary and it does not contribute to the 
Faddeev-Popov matrix JZ. Clearly, having a minimizing functional for the linear covariant gauge 
implies that the set of its local minima defines the first Gribov region £2 and that the corresponding 
Faddeev-Popov operator ^# is positive-definite. This should allow the extension of the Gribov- 
Zwanziger approach to the linear covariant gauge. In particular, one should be able to study the 
region Q. for the case of a gauge parameter E, ^ and to compare the results with the analytic study 



carried out in Reference [20] for small values of t,. 

In order to relate the lattice approach to the continuum, one should note [[if]] that the continuum 
relation (2T) can be made dimensionless by multiplying both sides by a 2 go. Since in d dimensions 
and in the SU(A^ C ) case one has j8 = 2N c /{a A ~ d g 2 ) ), it is clear that the lattice quantity 



P/(2N C 



a l g Q K b {x) 



x,b u 



1 

2^2 



I \a 2 g A b (x) 



x,b 



1 

2^2 



x,b 



(2.11) 



2 It is interesting to note that this functional can be interpreted as a spin-glass Hamiltonian for the spin variables g(x) 
with a random interaction given by (x) and with a random external magnetic field A(x). 
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gives 

^bV^x: = hl* x % [ A " w l 2 <2 - 12) 

in the fomial continuum limit. Thus, the continuum gauge parameter £ corresponds to a width 




(2.13) 



for the Gaussian distribution on the lattice and only for j8 = 2N C does one have a = y^. 



3. Numerical Simulations 



The functional Sucg\P 8 iSi^\ — see Eqs. (23) and (2JS) — is linear in the gauge transforma 



tions {g(x)}. Thus, the gauge-fixing algorithms usually employed in the Landau case [21, 22, 23] 
can be used also for the linear covariant gauge and, in principle, the numerical gauge fixing should 
not be a problem for this gauge. Nevertheless, as explained in the Abstract, any formulation of 
the linear covariant gauge on the lattice faces a conceptual problem. Indeed, the gluon field (x), 



usually defined as in Equation (2.6), is bounded while the functions A (x) satisfy a Gaussian dis- 



tribution [see Equation (2.2)] and are thus unbounded. Then, it is clear that Eq. (2.9) cannot be 
satisfied if A b (x) is too large [24]. As a consequence, one has to deal with convergence problems 
when numerically imposing the linear covariant gauge condition. Also, these problems become 
more severe as the width \J~E, of the Gaussian distribution becomes larger and/or as the lattice vol- 
ume becomes larger. On the lattice, as shown above, the width a of the Gaussian distribution is 
given by a = y / 2A r c ^/j3. Thus, these convergence problems are more severe also for small values 
of the coupling /3 . In particular, for j8 < 2N C the lattice width a is larger than the continuum width 

In References [[F7|, |||] we have presented tests of convergence of the numerical gauge fixing in 
the SU(2) case for relatively small lattice volumes, j8 = 4 and values of t, up to 0.5. [Recall that for 
j8 = 4 one has a = a/^ in the SU(2) case.] We have also checked that the quantity Di{p 2 )p 2 , where 
Dl(p 2 ) is the longitudinal gluon propagator, is approximately constant for all cases considered, as 
predicted by Slavnov-Taylor identities. This verification failed in previous formulations of the 
lattice linear covariant gauge [OL 1 6 1 . 

In order to overcome the convergence problems discussed above and be able to simulate at 
lattice coupling /3 smaller than 4 in the SU(2) case, we considered different discretizations of the 



gluon field. In particular, we used the angle (or logarithmic) projection [|26j and the stereographic 



projection [ |27[ ] (for a slightly different implementation of the stereographic projection see also 



[g8p). Note that, in the latter case, the gluon field is unbounded even for a finite lattice spacing a. 
Our results [ |7[ f?5j ] clearly show that the angle projection is already an improvement compared to 
the standard discretization and that the best convergence is obtained when using the stereographic 
projection. In References [|17L ^] we also presented preliminary results for the transverse gluon 
propagator using the stereographic projection. From these results one clearly sees that, as in Landau 
gauge, the transverse propagator is more infrared suppressed when the lattice volume increases. At 
the same time, for a fixed volume V, this propagator is also more infrared suppressed when the 
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gauge parameter £ increases. The latter result is in agreement with Reference |jT^]. One should, 
however, stress that the stereographic projection cannot be extended to SU(A/ C ) groups with N c > 2. 
Thus, in the SU(3) case one should probably rely on the logarithmic projection. 

Finally, one should note that, in the SU(2) case, the value a = \f%, i.e. j3 = 4, corresponds 
to a lattice spacing a as 0.001 fm. On the contrary, in the SU(3) case, one has a = \J~t, for j6 = 6, 
corresponding to a = 0.102 fm. Also, for a fixed t'Hooft coupling g^N c , we have j8 « N% and 
a oc y/ 1 fW c . This suggests that simulations for the linear covariant gauge are easier in the SU(iV c ) 
case for large N c . In Reference [|^] we tested this hypothesis by simulating the SU(2), SU(3) and 
SU(4) cases for a gauge parameter £ = 1 and lattice volumes up to V = 32 4 for several values of 
the coupling j3. We find that the convergence problems are indeed reduced when the number of 
colors N c is larger. 



4. The Limit £ ->■ 

In the continuum, the Landau gauge condition is defined by considering the usual Faddeev- 
Popov Lagrangian for the linear covariant gauge and by taking the limit E, — > 0. On the contrary, 
on the lattice, this gauge has been studied without considering this limit, but by imposing directly 
the gauge condition 

d^(x) = 0, (4.1) 
which is usually called the Lorenz gauge. The latter gauge condition is fixed numerically by min- 



imizing the functional ( |2.3| ). Using our implementation of the linear covariant gauge, it seems 
natural to study the Landau gauge by numerically considering the limit E, — > 0, in analogy with 
the definition in the continuum. One should also note that, in this limit, the width of the Gaussian 



distribution ( |2.2| ) goes to zero and, therefore, the convergence problems discussed above should be 
reduced — or eliminated — even for large lattice volumes and for j8 values in the scaling region. 
Moreover, since in the limit E, — > the Gaussian distribution becomes a Dirac delta function 5(A), 
this limit can be studied numerically by using different approximations of the delta function. 

Here we consider three possible sequences of functions, labelled by a parameter a, leading to 
a delta function in the limit a — > 0, i.e. 



1. the Gaussian distribution 

2. the Triangle distribution 

/r(A) = 

3. and the Rectangular distribution 

/«(A) = 



/g(A) 



e -A 2 /(2o£) 



(4.2) 



(l-A/«r) 
a T 

(l+A/« r ) 
a T 





for A G [0,a T ] 
forA G [-a T ,0] 
for A ^ [-a T ,a T ] 



(4.3) 



forA G [-Or/2, Or/2] 

forA^ [-a R /2,a R /2] 



(4.4) 



5 



Linear Covariant Gauge 



Attilio Cucchieri 



CM 



0.001 



0.00095 



0.0009 



0.00085 



0.0008 



0.00075 



0.0007 




0.5 1 



Figure 1: Plot of p 2 D L (p 2 ) as a function of p 2 (in lattice units) for a R = 0.1 (+), a T ~ 0.0707107 (□) and 
a G « 0.0288674 (A). We also show the theoretical value a 2 = aj/12 « 0.0008333. 



Note that the last two distributions are bounded. Also note that the average value of A 2 using these 
three distributions is respectively given by a 2 ,, and 0^/12. Thus, by setting 

V6 V12 

we have three different distributions with the same moment of inertia. 

We have done exploratory tests considering, in the SU(2) case, the lattice volume V = 32 4 
at j3 = 2.2. For the gauge-fixing parameters a we used the values GCg ~ 0.144338,0.0288674 
and 0.00288674 for the Gaussian distribution, a T « 0.353553,0.0707107 and 0.00707107 for the 
triangle distribution and a# = 0.5,0.1 and 0.01 for the rectangular distribution. Note that these 



values satisfy the relations in Eq. ( |4.5| ). Also, the values of Ofe correspond to the continuum values 



0.0114584,4.583310 4 and 4.583310 6 , respectively. For comparison, we have also done 



simulations directly at t, = 0, i.e. imposing the Lorenz condition (11). Results for the longitudinal 
Dl(p 2 ) and the transverse Dj(p 2 ) gluon propagators are shown in Figs. |]44[ One sees that the 
limit a — > is smoothly approached in the gluon sector and that the results are independent of the 
considered distribution. In particular, it is clear in Figs. that the theoretical prediction Di{p 2 ) = 
a 2 / p 2 , with a 2 = OCq = aj/12, is satisfied by the data obtained with these three distributions. At 
the same, a value of a 2 = CCq « 0. 144338 2 « 0.02 (see Fig. |j), which corresponds to the continuum 
value £ 0.01, already seems to give results in quantitative agreement with the Landau case. 
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Figure 2: Plot of Dl(p 2 ) as a function of p 2 (both in lattice units), using the logarithmic scale in both 
axes, for a R = 0.01 (+), a T « 0.00707107 (□) and a G « 0.00288674 (A). We also show the theoretical 
prediction a^/(12p 2 ). 
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Figure 3: Plot of Dj{p 2 ) as a function of p 2 (both in lattice units), using the logarithmic scale in both axes, 
for a R = 0.1 (+), a T ~ 0.0707107 (□) and a G w 0.0288674 (A). 
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16 p-i 1 r 




Figure 4: Plot of D T (p 2 ) as a function of p 2 (both in lattice units) for a G « 0.144338 (+), a G « 0.0288674 
(□), a G « 0.00288674 (A) and a G = (x). 

5. Conclusions 

We have found a minimizing functional for the linear covariant gauge that is a simple general- 
ization of the Landau-gauge functional. This approach solves most problems encountered in earlier 
implementations. Simulations for large lattice volumes, /3 values in the scaling region and a large 
gauge parameter £ can probably be done with SU(3) and SU(4) using the logarithmic projection, 
allowing a non-perturbative study of Green's functions. Finally, the approach to the limiting case 
£ — > 0, i.e. the Landau gauge, can also be studied numerically. Here we have investigated this limit 
considering three different distributions defining the gauge condition. 
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